A VARIATIONAL APPROACH TO DISLOCATION PROBLEMS 
FOR PERIODIC SCHRODINGER OPERATORS 



RAINER HEMPEL AND MARTIN KOHLMANN 

Abstract. As a simple model for lattice defects like grain boundaries in solid 
state physics we consider potentials which are obtained from a periodic poten- 
tial V = V(x, y) on R 2 with period lattice 1? by setting Wt(x, y) = V(x + t, y) 
for x < and Wt(x,y) = V(x,y) for x > 0, for t £ [0,1]. For Lipschitz- 
continuous V it is shown that the Schrodinger operators Ht = — A + Wt have 
spectrum (surface states) in the spectral gaps of Hq, for suitable t £ (0, 1). We 
also discuss the density of these surface states as compared to the density of 
the bulk. Our approach is variational and it is first applied to the well-known 
dislocation problem [Kl, K2] on the real line. We then proceed to the dis- 
location problem for an infinite strip and for the plane. In an appendix, we 
discuss regularity properties of the eigenvalue branches in the one-dimensional 
dislocation problem for suitable classes of potentials. 



1. Introduction 

In solid state physics, one first studies crystallized matter with a perfectly regular 
atomic structure where the atoms are located on a periodic lattice. However, most 
crystals are not perfectly periodic; in fact, the regular pattern of atoms may be 
disturbed by various defects which fall into two main classes: there are defects 
which leave the lattice unchanged (like impurities or vacancies), and there are 
more serious "geometric" defects of the lattice itself, cf. [AM], which may involve 
translations and rotation of portions of the lattice. Such lattice dislocations occur, 
in particular, at grain boundaries in alloys. These models are deterministic but 
may be generalized to include randomness. 

Many of the geometric defects mentioned above are accessible to mathematical 
analysis only after some idealization which leads to the following type of problem, 
cf. [DS]: there is a periodic potential V: R d — > R with period lattice 7L d and a 
Euclidean transformation T : R d — >• W l such that the potential coincides with V in 
the half-space {x £ R d | x\ > 0} and with V o T in {x\ < 0}. In the simplest cases 
T is translation in the direction of one of the coordinate axes, with again two main 
subcases: translation orthogonal to the hypcrplanc {x\ = 0} or translations that 
keep the a;i-coordinate fixed. In the present paper, we discuss the case d = 2 (where 
the coordinates are denoted by x and y) and we will mainly focus on translation 
in the a;-direction. In a forthcoming companion paper [HK] we will then study 
some aspects of the rotation problem where we take the given periodic potential 
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V in the right half-plane and a rotated version V o M§ in the left half-plane with 
M$ denoting rotation by the angle some results from the present paper will be 
essential for [HK]. 

The one-dimensional dislocation problem is particularly simple: Let V : R — > R 
be a periodic potential with period 1 and let 

for t e [0, 1]. The (self-adjoint) operator := + Wt is called the dislocation 

operator, t the dislocation parameter. There is quite a number of results available 
on this problem: it is well known and easy to see that the essential spectrum 
of H t does not depend on t for < t < 1; also H t cannot have any embedded 
eigenvalues. Furthermore, there is no singular continuous spectrum, cf. [DS]. For 
< t < 1, the operators H t may have bound states (discrete eigenvalues) located 
in the gaps of the essential spectrum. These eigenvalues and the corresponding 
resonances have been studied by Korotyaev [Kl, K2] in great detail, using powerful 
results from analytic function theory which are specific to the one-dimensional, 
periodic case. While, predictably, our results for the one-dimensional periodic case 
are weaker than Korotyaev's, our method of proof is very elementary and can be 
generalized in several directions; most importantly, we can apply our techniques to 
dislocation problems in dimensions greater than 1. In one dimension, we also give 
a more systematic treatment of regularity properties of the eigenvalue "branches" ; 
in particular, it is shown that the eigenvalue branches are Lipschitz-continuous if 

V is (locally) of bounded variation. 

The one-dimensional dislocation problem is mainly included to introduce and test 
our variational approach which is inspired by [DH, ADH] : we use approximations 
by problems on intervals (— n — t, n) with periodic boundary conditions where it is 
easy to control the spectral flow, and let n tend to 00. This idea can be adapted 
to the study of the translation problem for the strip E := R x (0, 1) in R 2 with 
periodic boundary conditions in the y- variable, say. In R 2 , we consider dislocation 
potentials W t defined by 

W t ( X ,y):=\lf' y Jt 1 "in (1-2) 

[ V(x + t, y), x < 0, 

for t G [0, 1]. On the strip S we obtain existence results for eigenvalues of St := 
— A + Wt in the spectral gaps of So- From that we easily derive that D t := — A + W t , 
acting in ^(K. 2 ), will have surface states with a non-zero density on an appropri- 
ate scale, for suitable t € (0, 1). To distinguish the bulk from the surface den- 
sity of states for this problem, we consider the operators —A + W t on squares 
Q n = (— n, n) 2 with Dirichlet boundary conditions, for n large, count the number 
of eigenvalues inside a compact subset of a non-degenerate spectral gap of Dq and 
scale with n~ 2 for the bulk and with n _1 for the surface states. Taking the limits 
n — > 00 (which exist as explained in [DS, EKSchrS]), we obtain the integrated den- 
sity of states measures pbuik(D t , I) for the bulk and p S urf(D t , J) for the surface states 
of this model; here I C K and J C R\o-(D ) are open intervals and J C R\a(D ). 
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Our main result can be described as follows: If (a, b) is a (non-trivial) spectral gap 
of the periodic operator —A + V, acting in ^(R 2 ), then for any compact interval 
[a, [i] C (a, b) with a < (3 there is a t € (0, 1) such that p sur f(D t , (a, (3)) > 0. Upper 
bounds for the surface density of states are discussed in [HK] . 

Our paper is organized as follows. Section 2 deals with dislocation on the real 
line. Here it is shown that the fc-th gap in the essential spectrum of H t (if it is 
open) is crossed by effectively k eigenvalues of H t as t increases from to 1. As an 
example, we discuss a periodic step potential in Section 3 where one can compute the 
eigenvalues of the dislocation operator numerically. Note that our calculations yield 
numbers which are exact up to finding the zeros of some transcendental functions. 
Related pictures can be found in [DPR] where a different numerical approach has 
been used. 

In Section 4 we adapt the method of Section 2 to the dislocation problem on the 
strip E. The results obtained for the strip then easily yield spectral information for 
the dislocation problem in the plane. Section 5 presents examples from the class of 
muffin tin potentials where one can "see" the motion of the eigenvalues rather 
directly for either translation in the x-direction or in the y-direction. Finally, 
we include an Appendix on regularity properties of the functions describing the 
eigenvalues of the dislocation operator H t in one dimension. 

For basic notation and definitions concerning self-adjoint operators in Hilbert 
space, we refer to [K, RS-I]. 

Acknowledgements. The authors would like to thank Andreas Ruschhaupt (Han- 
nover) and Evgeni Korotyaev (St. Petersburg) for helpful discussions. The authors 
are particularly indebted to Jiirgen Voigt, Dresden, who kindly contributed Lemma 
A. 6 in the Appendix. 



In this section, we study perturbations of periodic Schrddinger operators on the 
real line where the potential is obtained from a periodic potential by a coordinate 
shift on the left half-axis. 

Let ho denote the (unique) self-adjoint extension of — -^3 defined on C£°(R). 
Our basic class of potentials is given by 



Potentials V £ V belong to the class £i,ioc,unif(R) which coincides with the Kato- 
class on the real line; in particular, any V € V has relative form-bound zero with 
respect to h and thus the form-sum H of ho and V e V is well defined (cf. 
[CFrKS]). For V € V given, we define the dislocation potentials W t as in eqn. (1.1), 
for < t < 1; as before, the form-sum H t of ho and W t is well defined. 

We begin with some well-known results pertaining to the spectrum of H = Ho- 
As explained in [E, RS-IV], we have 



2. Dislocation on the real line 



V:={V £ L 1 ,i oc (R,R); Va; € M : V(x + 1) = V(x)} . 



(2.1) 



a(H)=a ess (H)=UZ 1 hk,l' k ], 



(2.2) 
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where the -fk and -f' k satisfy jk < "fl — 7fe+i> f° r au ^ <= N, and 7fc — > oo as k — > oo. 
Moreover, the spectrum of is purely absolutely continuous. The intervals [7^, 7^] 
are called the spectral bands of H. The open intervals Tk := (7k)7fc+i) are the 
spectral gaps of if; we say the fc-th gap is open or non-degenerate if 7^+1 > 7^,. 

In order to determine the essential spectrum of H t for < t < 1, we introduce 
Dirichlet boundary conditions at x = for the operator H and at x — and 
x = — t for H t to obtain the operators 

H D = H-(BH+, H t!D = H t -®H { _ tfi) ®H+, (2.3) 

where H ± acts in with a Dirichlet boundary condition at 0, in (—00, —t) 
with Dirichlet boundary condition at —t and Ht t ^ in (— i, 0) with Dirichlet 
boundary conditions at —t and 0. Since iI(_ ti o) has purely discrete spectrum 
and since the operators and if - are unitarily equivalent, we conclude that 
aessiHo) = cr ess (£/*,£>)■ It is well known that decoupling by (a finite number of) 
Dirichlet boundary conditions leads to compact perturbations of the corresponding 
resolvents (in fact, perturbations of finite rank) and thus Weyl's essential spectrum 
theorem yields cr css (H D ) = cr css (H) and <7 css (H t , D ) = a ess (H t ). 

In addition to the essential spectrum, the operators H t may have discrete eigen- 
values below the infimum of the essential spectrum and inside any (non-degenerate) 
gap, for t € (0, 1); these eigenvalues are simple. The eigenvalues of H t in the gaps 
of H depend continuously on t; cf. the Appendix for a brief exposition of the rel- 
evant perturbational arguments, which are fairly standard. A more complete and 
precise picture is established in the following lemma which says that the discrete 
eigenvalues of Ht inside a given gap Tk of H can be described by an (at most) count- 
able, locally finite family of continuous functions, defined on suitable subintervals 
of [0,1]. 

Lemma 2.1. Let k <G N and suppose that the gap Tk of H is open, i.e., 7^ < 7fe+i- 
Then there is a (finite or countable) family of continuous functions fj : (aj,f3j) — > 
Tk, where < aj < f3j < 1, with the following properties: 

(i) fj(t) is an eigenvalue of H t , for all aj < t < f3j and for all j. Conversely, 
for any t G (0, 1) and any eigenvalue E e Tj. of H t there is a unique index 
j such that fj(t) = E. 

(ii) As t I aj (or t t Pj)> the limit of fj(t) exists and belongs to the set 
WkHk+i}- 

(Hi) For all but a finite number of indices j the range of fj does not intersect a 
given compact subinterval [a',b'] C IV 

For the convenience of the reader, we include a proof in the Appendix. Under 
stronger assumptions on V one can show that the eigenvalue branches are Holder- 
or Lipschitz-continuous, or even analytic; cf. the Appendix. Additional information 
on the eigenvalue functions fj can be found in [Kl, K2]. 

It is our aim in this section to show that at least k eigenvalues move from 
the upper to the lower edge of the fc-th gap as the dislocation parameter ranges 
from to 1. Using the notation of Lemma 2.1 and writing /i(a,) := lim t _|_ Qi fi(t), 
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fi((3i) := limtt^ fi(t), wc now define 

N k := #{i | Ma,) = 7fe+1 , /.(ft) = 7^} - #{i | / t (a s ) = 7^, = 7fe+1 }. 

(2.4) 

Thus A/fe is precisely the number of eigenvalue branches of H t that cross the k- 
th gap moving from the upper to the lower edge minus the number crossing from 
the lower to the upper edge. Put differently, A4 is the spectral multiplicity which 
effectively crosses the gap Tk in downwards direction as t increases from to 1. 
Our main result in this section says that A4 = k, provided the fc-th gap is open: 

Theorem 2.2. Let VeP and suppose that the k-th spectral gap of H is open, i.e., 
Ik < 7fe+i- ThenAfk = k. 

Again, the results obtained by Korotyaev in [Kl, K2] are more detailed; e.g., it 
is shown that, for any t £ (0, 1), the dislocation operator H t has two unique states 
(an eigenvalue and a resonance) in any given gap of the periodic problem. On the 
other hand, our variational arguments are more flexible and allow an extension to 
higher dimensions, as will be seen in the sequel. In this sense, the importance of 
this section lies in testing our approach in the simplest possible case. For further 
reading concerning the spectral flow through the gaps of perturbed Schrodinger 
operators, we recommend [P, Saf]. 

The main idea of our proof — somewhat reminiscent of [DH, ADH] — goes as fol- 
lows: consider a sequence of approximations on intervals (— n — t, n) with associated 
,2 

operators H n .t — — ^3 + W 7 * with periodic boundary conditions. We first observe 
that the gap Tk is free of eigenvalues of H n $ and H n ^ since both operators are ob- 
tained by restricting a periodic operator on the real line to some interval of length 
equal to an entire multiple of the period, with periodic boundary conditions. Sec- 
ond, the operators H n j have purely discrete spectrum and it follows from Floquet 
theory (cf. [E, RS-IV]) that H n $ has precisely 2n eigenvalues in each band while 
H n \ has precisely In + 1 eigenvalues in each band. As a consequence, effectively 
k eigenvalues of £f„ jt must cross any fixed E e Tk as t goes from to 1. To obtain 
the result of Theorem 2.2 we only have to take the limit n — > 00. Here we employ 
several technical lemmas. In the first one, we show that the eigenvalues of the 
family H n , t depend continuously on the dislocation parameter. 

Lemma 2.3. The eigenvalues of if„ it depend continuously on t e [0, 1]. 

Proof. We may assume that the eigenvalues of H n ^ are numbered according to 
min-max. Since the Hilbert space L^—n — t, n) depends on t, we use the unitary 
mappings 

U n ,t- L 2 (-n-t,n) L 2 (-n,n), (U n , t f)(x) := ^/b\^f(a n ^ t x), (2.5) 

where a„ it := Let H n ^ t := U n . t H n jU~j and W t (x) := W t (a ntt x) so that 

(writing a = a nyt ) 

H n ,t = °~~ 2 ho + W t {x) = a- 2 (h + <J 2 W t {x)). (2.6) 

It is easy to see that the mapping [0, 1] 3 (4 cr 2 W t € Li{—n,n) is continuous. 
Now the usual perturbational and variational arguments for quadratic forms ([K] 
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and the Appendix) imply that the eigenvalues of ho + <? 2 Wt depend continuously 
on t, and then the same is true for the eigenvalues of H njt . □ 

The next lemma is to establish a connection between the spectra of H t and H n t 
for < t < 1 and n large. In the proof and henceforth, we will make use of the 
following cut-off functions: We pick some ip e C^°(— 2,2) with < ip < 1 and 
<p(x) = 1 for \x\ < 1. For k e (0, oo) we then define (fk{x) '■— <p{x/k) so that 
supp (p k C {-2k, 2k), (p k (x) = 1 for |a;| < k, \(p' k (x)\ < Ck' 1 and \<p£(x)\ < Ck~ 2 . 
Finally, we let ^ := 1 — ipk- F° r an y self-adjoint operator T we denote the spectral 
projection associated with an interval IcRby Pi(T) and we write dim Pi{T) to 
denote the dimension of the range of the projection Pi(T). 

Lemma 2.4. Let fceN with Tk ^ 0. Let t e (0, 1) and suppose that r\\ < ry 2 e Tk 
are such that r\\,r\i £ a(H t ). Then there is an n e N such that 771,7/2 a(H nit ) 
for n> n , and 

dim P( Vum) (H t ) = dim P (mim) (H ntt ), n > n . (2.7) 

Proof. In the subsequent calculations, we always take k :— ra/4, for nefl. 

(1) Let E E (T71 , 772) n o-(Ht) with associated normalized eigenfunction u. Then 
Uk ■= ifiku E D(H n .t), Hn.tUk — H t Uk and \uk\ — > 1 as n — > 00. Since 

\\H n . t u k -Eu k \\<2- \\u>\\ + l^'IU \\u\\ , (2.8) 

it is now easy to conclude that dim P( Vim )(H nit ) > dim P( m ,r) 2 )(Ht) for n large. 

(2) We next assume for a contradiction that 77 e Tk satisfies 77 € cr(H n j) for 
infinitely many neff. Then there is a subsequence (rij)j e N C N s.th. 77 e a(H njit ); 
we let u nj! i G D(H njtt ) denote a normalized eigenfunction and set 

Vl, nj ■ = tPkjUnj.U V 2 . n] ■= IpkjUn^t, (2.9) 

so that v\. nj € D(H t ) and \\(H t — v) v i,n j \\ — > as j — > oo by a similar estimate as 
in part (1) (and using a simple bound for t || which follows from the fact that 
V has relative form-bound zero w.r.t. ho.) Let us now show that V2, nj — > (and 
hence ||ui )TlJ || — > 1) as j — > oo: The function 

belongs to the domain of ff nj ,o and H n]fi V2, n] = [H nj>t V2, nj ]~ , where [-]~ is defined 
in analogy with eqn. (2.10). Since we also have (H njtt — r))v2, nj — > 0, as j — » oo, 
we see that (H nj , — 77)^2, «j — > 0. But dist(?7, a(H n fi)) > S > for all n and the 
Spectral Theorem implies that ||u2,n^ || ~~ ^ as j — > oo. We have thus shown that 
«j || - ^ 1 an d \(Ht — fi)vi, nj || — > which implies that 77 e cr(H t ). 

(3) It remains to show that dim P^ i m ^{H n ^ t ) < dim Pt m t7]2 \(H t ) , for n large. 
The proof by contradiction follows the lines of part (2); instead of a sequence of 
functions u nj we work with an orthonormal system Un] , . . . , Un] of eigenfunctions 
where I = dim P( m , m ){Ht + !)• We leave the details to the reader. □ 
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Remark 2.5. In fact, using standard exponential decay estimates for resolvents of 
Schrodinger operators, cf. [S], it can be shown that the eigenvalues of H t and H n t 
in the gap Tk are exponentially close, for n large; e.g., if E S cr(H t ) n Tk for some 
t G (0, 1), then there are constants c > and a > s.th. the operators H n>t have 
an eigenvalue in (E — ce~ an , E + cc~ an ), for n large. There is a similar converse 
statement with the roles of H t and H n t exchanged; cf. also Remark 4.2 for further 
discussion. 

The desired connection between the spectral flow for (-ff n ,t)o<t<i an d (H t )o<t<i 
is obtained by applying Lemma 2.4 at suitable U € [0, 1] and ryi,, < r\2,i € IV Wc 
now construct an appropriate partition of the parameter interval [0, 1]. 

Lemma 2.6. Let fceN with T k ^ 0. Then there exists a partition = t < t\ < 
■ ■ ■ < tx-i <tK = 1 and there exist Ej G Tk and n G N such that 

Ej i <j(H t )U<j(H n . t ), Vt £ [tj-utj], j = l,...,K, n>n . (2.11) 

Proof. For any t € [0,1] there exists r\ t € Tfe such that % ^ o-(H t ). Since the 
spectrum of £f t depends continuously on the parameter t there also exists e — e t > 
such that 774 ^ cr(H T ) for all t e (t — St,t + et). By compactness, we can find a 
partition {Tj)o<j<K (with to = 0, tk = 1) such that the intervals (t^ — £j,Tj + £j) 
cover [0, 1]. Set Ej := r\ Ti . We next pick arbitrary points tj € (rj, Tj + £,•)□ {tj + i — 
£j+i,Tj+i), for j = 1, . . . ,K — 1, set i = 0, t# = 1 and see that Ej £ cr(H t ) for 
tj-i < t < tj, j = 1,...,K. By Lemma 2.4, using Lemma 2.3 combined with a 
simple compactness argument, we then find that we also have Ej ^ o~(H n . t ) for 
t G [tj-i,tj] and n large. □ 

We are now ready for the proof of Theorem 2.2. 

Proof of Theorem 2.2. Let Ej be as in Lemma 2.6 and A4 as in eqn. (2.4). We fix 
some E e T k such that E > Ej for j = 0, . . . , K and E <£ a(H tj ) U o-{H n . tj ) for 
j = 0, . . . , K and for all n large. It is then easy to see that 

K 

Afk = (dim P {Ej ^(H tj ) - dim P {Ej ft{H tj _S) (2.12) 
j'=i 

and that 

dim P (-oo„E)(#«,i) - dim P {-ooM Hn ^ 

= J2 ( dim PiE^E^Hnti) - dim P (Bj> B)(^n,t J _i)^- 

J=1 (2.13) 
The LHS of (2.13) equation equals fc. Furthermore, by Lemma 2.4, we have 

dim P ( Ej M H ^ = dim P (£„s)(^) (2- 14 ) 
for all j and all n large, and the desired result follows. □ 
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3. A ONE-DIMENSIONAL PERIODIC STEP POTENTIAL 

In this section, we study the one-dimensional 27r-periodic potential 

(While the other sections of this paper deal with 1-periodic potentials, we have 
preferred to work here with period 2ir in order to keep the explicit calculations done 
by hand as simple as possible.) To obtain the band-gap structure of H = —-^z + V, 
we compute the discriminant function 

where tpi(-;E) and Lp 2 (- ;E) solve the equation 

-it" + (V - E)u = (3.3) 

and satisfy the boundary conditions 

ipi(0;E) = ip' 2 (0;E) = 1 and ^(0; E) = ip 2 (0; E) = 0. (3.4) 

The matrix M{E) on the RHS of (3.2) is called the monodromy matrix. A simple 
computation shows that [—1/2, 1/2] C Ti, where Ti is the first spectral gap of H 
(with numbering according to Floquet theory). Note that the gap edges of I\ also 
equal the first eigenvalue in the (semi-)periodic eigenvalue problem for —3^2 + V 
in 1,2(0,271"), cf., e.g., [E, CL]. 

As explained in [E, RS-IV], for any E ^ cr(H), there are two solutions ip±(x; E) e 
C 1 (R), square integrable at ±00, of (3.3); in fact, the functions (p±(x;E) are ex- 
ponentially decaying at ±00 and exponentially increasing at =Foo. In our example, 
the dislocation potential Wt for t € (0, 1) will produce a bound state at E if and 
only if the boundary conditions coming from f+(0; E) and <p-(t; E) match up, i.e., 



An equivalent condition for (3.5) is the equality of the ratio functions ^, and 

y+(Q;£) 



<P-(t;E) = <p+(0;E) and <p'_(t; E) = ip' + (0; E). (3.5) 

f-(t;E) 
•p'-(t;E) 

y^"(0;gj ' PPR]- We compute the Floquet solutions ip± by solving the equation 
-u" + (V - E)u = for x < and x > and for E varying in [-1/2,1/2], 
assuming that (u(0), w'(0)) equals an appropriate eigenvector of M(E). Note that, 
since D(E) < —2, both eigenvalues of M(E) are negative and not equal to — 1. 
Finally, we divide [—1/2,1/2] into 100 subintervals of equal length and compute 
numerical values for t such that 



V -{t;E) ^+(0;E) 



<p'_(t;E) <p' + (0;E) 



< e, (3.6) 



where the error e > is small. This leads to the following plot of t M> E(t), see 
Fig. 1. 
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Figure 1: An eigenvalue branch of H t in the first spectral gap. 



4. Periodic potentials on the strip and the plane 



Let V: M 2 — > M be Z 2 -periodic and Lipschitz-continuous and let £ = M x (0, 1) 
denote the infinite strip of width 1. We denote by S t the (self-adjoint) operator 
—A + Wt, acting in L2CS), with periodic boundary conditions in the y-variable and 
with Wt defined as in cqn. (1.2); again, the parameter t ranges between and 1. 
Since So is periodic in the x-variablc, its spectrum has a band-gap structure. 

We first observe that the essential spectrum of the family St does not depend 
on the parameter t, i.e., a ess (St) — ct css (Sq) for all t G [0,1]. As in Section 2, 
this follows from the compactness of (St — c) _1 — (S t .D — c )~\ where St,D is St 
with an additional Dirichlet boundary condition at x = 0, say. (While, in one 
dimension, adding in a Dirichlet boundary condition at a single point causes a rank- 
one perturbation of the resolvent, the resolvent difference is now Hilbcrt- Schmidt, 
which can be seen from the following well-known line of argument: If — As denotes 
the (negative) Laplacian in L 2 (T,) and — As ; d is the (negative) Laplacian in L 2 (T<) 
with an additional Dirichlet boundary condition at x = 0, then (—As + 1) _1 — 
(— Ay, ; d + l) -1 has an integral kernel which can be written down explicitly using 
the Green's function for —As and the reflection principle.) 

While the essential spectrum of the family St does not change as t ranges through 
[0, 1], St will have discrete eigenvalues in the spectral gaps of So for appropriate 
values of t. We have the following result. 

Theorem 4.1. Let (a,b), a < b, denote a spectral gap of St and let E G (a, b). 
Then there exists t = tE € (0, 1) such that E is a discrete eigenvalue of St- 

Proof. (1) As on the real line, we work with approximating problems on finite size 
sections of the infinite strip S. Let 



and consider S n . t := — A + W t acting in L 2 (S„ it ) with periodic boundary conditions 
in both coordinates. The operator S n ^ has compact resolvent and purely discrete 
spectrum accumulating only at +00. The rectangles S„ i0 (respectively, S n ,i) consist 



S„, t := (-n-t,n) x (0,1), n e N, 



(4.1) 
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of 2n (respectively, 2n + 1) period cells. By routine arguments (see, e.g., [RS-IV, 
E]), the number of eigenvalues below the gap (a, b) is an integer multiple of the 
number of cells in these rectangles; we conclude, that eigenvalues of S n j must cross 
the gap as t increases from to 1. 

(2) Let E G (a, b). According to (1), for any n G N we can find t n G (0, 1) such that 
E G o _ disc('S'n,t„)i then there are eigenfunctions u n G D(S , „ ;trl ) with S n j n u n = Eu n , 
\\u n \\ = 1, and ||Vu„| < C for some constant C > 0. We now choose cut-off functions 
(p n as in Section 2 and denote the natural extension to M 2 again by <p n . We also 
let ip n = 1 — ip n . Clearly, 

)||, (S n , tB -£7)(V )\\<c/n, (4.2) 

for some c > 0. There is a subsequence [t nj )jeN C (i n ) n eN and t G [0,1] s.th. 
t nj -^t&sj^r oo. Since is Lipschitz, we may infer from (4.2) that 

rij /A^Uj 

) \\->0, i^oo, (4.3) 

and it remains to show that |^ n /4U„|| — > so that ||</9„/4U„|| — > 1. We associate 
with functions w : E nj t — > C functions £ : X„,o — > C by 

in analogy with eqn. (2.10). Then [VV4 U «]~ € D(S n> o) and 

||(<SVi,o - E)[ip n / 4 u n }~\\ = \\{S n . tn - E)(ip n / 4 u n )\\ < c/n. (4.5) 

Since (a, b) (1 cr(S n ,o) = f° r au ^ G N, and since £7 G (a, 6), the Spectral Theorem 
implies that [ip n /4U n ]~ — > (and therefore also ip n /iU n — > 0) as n — > oo. 

We therefore have shown that the functions ti nj . := <p n /iU nj for j G N satisfy 
I (Sj — E)v nj || —¥ and || — >• 1 as j — > oo which implies E G o-(Sj). □ 

Remark 4.2. By a well-known line of argument, one can obtain exponential lo- 
calization of the eigenfunctions of St near the interface {(x, y) \ x = 0}. Since we 
will use exponential localization in a more systematic way in the forthcoming paper 
[HK] we only give a brief sketch here: Suppose that E G (a, b) and t G (0, 1) satisfy 
£7 G cr(St). Let it G D(So) = D(S t ) denote a normalized cigenfunction and let ip n , 
n G N, be as in the proof of Theorem 4.1. As above, we have 

(St - E)(tp n u) = -2Vip n ■ Vu - (A<p n )u =: r n , (4.6) 

where |r„| < c/n, for n G N. Since r n has support in the interval (— 2n — 1, 2n) we 
now see that there exist constants C > and a > such that 

||XM>4„HI < ||x|x|>4n(S , t - E)- 1 ^]] < Ce~ an , (4.7) 

by standard exponential decay estimates for the resolvent kernel of Schrodingcr 
operators (cf., e.g., [S], [HK]). 

We now turn to the dislocation problem on the plane R 2 where we study the 
operators 

D t = -A + W t , 0<t<l. (4.8) 
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Denote by St($) the operator St with ^-periodic boundary conditions in the y- 
variable. Since W t is periodic with respect to y, we have 




and hence the spectrum of D t has a band-gap structure; furthermore, D t has no 
singular continuous part, cf. [DS, FS]. As for the spectrum of S t inside the gaps of 
So, Theorem 4.1 leads to the following result. 

Theorem 4.3. Let (a,b) denote a spectral gap of Do, a > inf a css (Do), and let 
E € (a, b). Then there exists t = tE € (0, 1) with E G a(D t ). 

Proof. Let ip n u n G D(S t ) as in part (2) of the proof of Theorem 4.1 denote an ap- 
proximate solution of the eigenvalue problem for St and E. We extend u n to a func- 
tion u n (x,y) on R 2 which is periodic in y. Writing $„ = & n (x,y) := ip n (x)^ n (y) 
we compute 

(D t - E)(<5> n u n ) = (-d 2 x -d 2 y + W t -E) (<p n (x)y n (y)u n (x,y)) 

= <Pn{y)[(S t - E)(<p n (x)u n )]~ - ipn(x) (2ip' n {y)d y u n + ip'n(y)u n ) . 

(4.10) 

The norms of the three terms on the RHS can be estimated (up to a constant which 
is independent of n) by en, and ^n, respectively, and we see that 

\\(D t - E)(^ n u n )l < c (l + ne), (4.11) 

while |$„M n | > con with a constant Co > 0. This implies the desired result. □ 

Remark 4.4. We learn from the above proof that there are functions 

v n = v n (x, y) := * A > n u n (4.12) 

that satisfy \\v n \\ = 1, supp v n C [— n,n] 2 and 

(D t - E)v n ->• 0, n -> oo. (4.13) 

These functions play a key role in our analysis of the rotation problem at small 
angle in [HK]. 

We finally turn to a brief discussion of the i.d.s. (the integrated density of states 
[V]) for the dislocation operators D t . We adopt the natural distinction of [DS, 
EKSchrS, KS] between bulk and surface states. Roughly speaking, the bulk states 
correspond to states away from the interface with energies in the spectral bands 
while the surface states for < t < 1 are produced by the interface and are (expo- 
nentially) localized near the interface. The (integrated) density of states measures 
for the bulk and surface states use a different scaling factor in the following defini- 
tion: restricting D t to large squares Q n = (—n,n) 2 and taking Dirichlet boundary 
conditions, we obtain the operators D^ n \ For / clan open interval, let N(I, D^) 
denote the number of eigenvalues of D^ in /, counting multiplicities. We then de- 
fine for open intervals / cR and J C R \ cr(D ) with J C M \ a(D ) 

Pbuik(/,A)= lim -^N(I,D[ n) ), p surf (J,A)= lim ^N( J, D t (n) ). (4.14) 
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The existence of the limits in (4.14) has been established in [EKSchrS, KS] for 
crgodic Schrodinger operators. Note that the surface density of states measure is 
defined (and possibly non-zero) for subintervals of the spectral bands, but then 
cqn. (4.14) is not suited to capture the surface states (cf. [EKSchrS, KS]). 

The fact that the surface density of states exists does not mean it is non-zero 
and there are only rare examples where we know p sur f to be non-trivial. It is one 
of the main results of the present paper to show that dislocation moves enough 
states through the gap to have a non-trivial surface density of states, for suitable 
parameters t. Indeed, it is now easy to derive the following result: 

Corollary 4.5. Let (a, b) be a spectral gap of Dq with a > inf er css (_Do), and let ^ 
J C (a, b) be an open interval. Then there is at e (0, 1) such that p sur f(J, Dt) > 0. 

Proof. Let [a, (3} C J with a</5,fix£e (a,/3), and let < e < mm{E-a, /3-E}. 
By Theorem 4.3 and Remark 4.4 there exist t = ts € (0, 1) and a function uq in 
the domain of D t satisfying ||ito|| = 1, supp uq compact, and \\(D t — E)uo\ < e. 
Let v e N be such that supp w C (— v, v) 2 \ note that, in the present proof, v 
corresponds to the n of Remark 4.4. We then see that the functions cpt, defined by 
tpk(x, y) := uq(x, y — 2ku) for k € N, have pairwise disjoint supports, are of norm 1, 
and satisfy \\(D t — E)ipk\ < £■ Furthermore, we have supp ipk C (— n,n) 2 provided 
(2k + l)v < n. Denoting M n '■= span-j^ | k e N, k < — 1)}, it is clear 
that dim.Mn > n/(3z^), for all n large. Let Af n denote the range of the spectral 
projection P( Q ,. / g)(I?("' ) ) of associated with the interval (a,/3); we will show that 
dimA/"„ > dim A4 n which implies the desired result. If we assume for a contradiction 
that dimA/"„ < dim A4 n for some n € N, we can find a function v £ A4 n fl N,^ of 
norm 1. By the Spectral Theorem, ||(£*j™^ — E)v\ > e. On the other hand, v is a 
finite linear combination of the tpk, which implies |(-Di™' ) — E)v\\ < e. □ 

We will continue the discussion of bulk versus surface states in the companion 
paper [HK] where a corresponding upper bound of the form N(J, D^) < cn log n 
is provided. 

5. Muffin tin potentials 

Here we present some simple examples where one can see the behavior of surface 
states directly. We will deal with Z 2 -periodic muffin tin potentials of infinite height 
(or depth) on the plane R 2 which can be specified by fixing a radius < r < 1/2 
for the discs where the potential vanishes, and the center P — (x , yo) € [0, l) 2 for 
the generic disc. In other words, we consider the periodic sets 

£l r ,p := U( ii _ 7 -) eZ 2B r (Po + 

and we let V = V r .p a be zero on f2 r .p while we assume that V is infinite on 
R 2 \ Q r ,p - If Hij is the Dirichlet Laplacian of the disc B r (Po + then the 

form-sum of —A and V r ,p is ®(i.j) e z 2 Hij ■ Without loss of generality, we may as- 
sume yo = henceforth. 
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(1) Dislocation in the x-direction. Here muffin tin potentials yield an illus- 
tration for some of the phenomena encountered in Section 4. In the simplest case 
we would take x = 1/2 so that the disks B r (l/2 + for i e No and j e Z, 
will not intersect or touch the interface {(x,y) | x — 0}. Defining the dislocation 
potential W t as in Section 4, we see that there are bulk states given by the Dirichlct 
eigenvalues of all the discs that do not meet the interface, and there may be surface 
states given as the Dirichlet eigenvalues of the sets B r (l/2 — t,j) n {x < 0} for 
j'eZandl/2-r<t<l/2 + r. 

More precisely, let fik — Hk(r) denote the Dirichlet eigenvalues of the Laplacian 
on the disc of radius r, ordered by min-max and repeated according to their respec- 
tive multiplicities. The Dirichlet eigenvalues of the domains B r (l/2 — t, 0)n{x < 0}, 
1/2 — r < t < 1/2 + r, are denoted as Afc(t) = Xk{t, r); they are continuous, mono- 
tonically decreasing functions of t and converge to Hk as t j" 1/2 + r and to +oo as 
t 4- 1/2 — r. In this simple model, the eigenvalues correspond to the bands of a 
periodic operator. We see that the gaps are crossed by surface states as t increases 
from to 1, in accordance with the results of Section 4 (Corollary 4.5). 

Along the same lines, one can easily analyze examples where xq is different from 
1/2; here more complicated geometric shapes may come into play. In [HK] we will 
again use muffin tin potentials as examples for the rotation problem. In that paper, 
we will also discuss approximations by muffin tin potentials of height n and their 
limit as n — > oo. 

(2) Dislocation in the y-direction. This problem has not been considered so 
far. We include a brief discussion of this case for two reasons: on the one side, we 
observe a new phenomenon which did not appear so far; on the other hand, one can 
see from our example that, presumably, there is no general theorem for translation 
of the left half-plane in the y-direction. 

Let V — V r denote the muffin tin potential defined above, with x = yo = 0. 
We then let W t coincide with V in the right half-plane, while we take W t (x, y) — 
V(x, y — t) in the left half-plane. At the interface {x = 0} we see half-discs on the 
left and on the right with the half-discs on the right being fixed while the half-discs 
on the left are shifted by t in the y-direction. The surface states correspond to the 
states of the Dirichlet Laplacian on the union O tjr;sur f of these half-discs. There 
are two cases: either f2( jr;S urf is connected and we have a scattering channel along 
the interface, or rit ir;sur f is the disjoint union of a sequence of bounded domains; cf. 
Figure 2. In the second case, the eigenvalues on such domains start at the Dirichlet 
eigenvalues of the disc of radius r, increase up to the corresponding eigenvalues of 
a half-disc, and then move down again to where they started. For 1/4 < r < 1/2, 
the picture is more complicated: If we let t = 1 — 2r, n = 2r, we find that the sets 
^t,r;surf are disconnected for < t < To and for n < t < 1; for To < t < n, however, 
^t,r;surf is connected and forms a periodic wave guide with purely a.c. spectrum 
[SW] ; cf. also [DS] . We therefore observe a dramatic change in the spectrum of the 
dislocation operators: for t G [0, To] U [n, 1] the surface states in the gap are given 
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by eigenvalues of infinite multiplicity while for t € (to,ti) the surface states form 
bands of a.c. spectrum in the gaps. 




Figure 2: Muffin tins: two cases for dislocation in the y-direction. 



Note that, if we had chosen xo = 1/2, then nothing at all would have happened 
for translation in the y-direction. 

(The authors thank A. Ruschhaupt, Hannover, for asking about translation in 
the y-direction.) 

6. Appendix: continuity and regularity of eigenvalues 

In this appendix, we discuss several basic facts concerning continuity and regu- 
larity of the eigenvalue branches for the one-dimensional dislocation problem. We 
hrst consider potentials V from the class V C Li,ioc,unif(R) as in (2.1) where the 
eigenvalues are continuous functions of the dislocation parameter t. In the subse- 
quent estimates we will use 

l^ll.locunif^^p f V+1 \V(x)\dx (Al) 

as a natural norm on Li,ioc,unif(R). As is well known (cf., e.g., [CFrKS]), any 
potential V G Li.ioc,unif(BO is relatively form-bounded with respect to ho with 
relative form-bound zero. More precisely, we have the following lemma. 

A.l. Lemma. For any e > there exists a constant C e > such that for any 
V € £i,ioc,unifW we have 

/|V|M 2 d*< ^U.^ (e ||^f + C £ |M| 2 ) , ^tf(l). (A2) 

Proof. For / G C^°(R) with support contained in (0,e) we have l/l^ < v^ll/'ll- 
Let (Cn)neN denote a (locally finite) partition of unity on the real line with the 
properties: supp Ci C (0, e), each is a translate of Ci, M := sup.j, gR S«eN ICOf)) 2 
is finite and J2 n eN^n( x ) — 1 f° r au ^ G K. By the IMS localization formula (see 
[CFrKS]), we have for any ip G C C °°(M), 

OO OO OO 

ii^f = (-<p», V ) = £ KCn^f - E ion 2 > E KCn^f - M imi 2 . 

n— 1 ri—1 n — 1 



R. HEMPEL, M. KOHLMANN 



15 



so that 

in«)ii^)i 2 d«<Ei^i~ / w»oid* 

„=1 "'supp Cn 

l,loc,unif * 

The general case follows by approximation and Fatou's lemma. □ 
For V e V, the function 

ti v (s):=[ \V(x + s)-V(x)\dx, < s < 1, (A3) 
Jo 

is continuous and iV(s) — > 0, as s — > 0. Furthermore, for W t is as in eqn. (1.1), we 
have 1 W t — W t > [ x loc unif = #y(i — £')■ This leads to the following lemma. 

A. 2. Lemma. Let V e V, E Q e R \ cr{H to ), and write e := dist(E , a(H to )). 
Then there is r > such that H t has no spectrum in (E — £o/2, E + £o/2) for 
\t — to | < T o- Furthermore, there exists a constant C > smc/i t/iat /or some 

Tl G (0,T ) 

\\iHt-Eo)- 1 -{H to -E )- l \\<m v {t-t ), \t-t \<n. (AA) 

Proof. Without loss of generality we may assume that V > 1. Let h t denote the 
quadratic form associated with H t . Applying Lemma A.l (with e := 1) we see that 

KM-htMl < f iWt-w^WufdxKddvit-to^M u€H\m), 

with some constant C\. The desired result now follows by [K; Thm. VI-3.9]. □ 

We therefore see that H tn — > H to in the sense of norm resolvent convergence if 
to € [0, 1], (t n )nem C [0, 1] and t n — > t . By standard arguments, this implies that 
the discrete eigenvalues of H t depend continuously on t. We are now prepared for 
the proof of Lemma 2.1. 

Proof of Lemma 2 A. We consider t € T, the flat one-dimensional torus, and we 
denote the spectral gap by (a, 6). Let [a', 6'] C (a, 6). 

(1) Let (77, t) e (a, 6) x T. Since cr(H T ) n (a, 6) is a discrete set, and since cr(H t ) 
depends continuously on t, there is a neighborhood U V , T C (a, 6) x T of (77, r) of the 
form L/^t = (771,772) x (ti,t 2 ) belonging to either of the two following types: 

Type (1): For t\ < t < r 2 we have a(iJ t ) n (771,772) = 0, 

or 

Type (2): 77 is an eigenvalue of H T and there is a continuous function / : (ti,t 2 ) — > 
{V11II2) suc h that /(t) is an eigenvalue of -ff t ; iJ t has no further eigenvalues in 
(»7i,%), for ti < f < r 2 . 

Now the family {t/^r | (77, r) G (a, 6) x T} is an open cover of (a, &) x T and there 
exists a finite selection {J7 w , ri }i=i,...,jv such that 

[a', 6'] xTcuf^ j)Ti . 
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As a first consequence, we see that there is at most a finite number of functions 
that describe the spectrum of H t in the open set uf^l^.^ D [a', b'] x T. 

(2) Suppose that (77, r) e (a, b) x T is such that rj € <r(H T ) and let /: (ti,t 2 ) — > 
(771, 1J2) as above. Consider a sequence (tj)jeN C (n, r 2 ) with tj — > n. We can find 
a subsequence (tj fc )/ceN such that f{tj k ) — > fj for some 77 G , ^2] - It is easy to see 
that fj € a(H T1 ). If f\ e (a, 6) the point (77, Ti) has a neighborhood L^ jT1 of type (2) 
and we can extend the domain of definition of / beyond n . It follows that there 
exist a maximal open interval (a,/3) C (0,1) and a (unique) continuous extension 
/: (a, f3) — > (a, 6) of / such that f(t) is an eigenvalue of Ht for all t e (a, /3). 

(3) It remains to show that f(t) converges to a band edge as t 1 ct and as rj t By 
the same argument as above, we find that any sequence (ij)jgN C (a, /3) satisfying 
tj — > a has a subsequence (tj k )k&>s such that f(tj k ) — > 77 for some 77 G [a, &]. Here 
77 ^ (a, 6) because otherwise we could again extend the domain of definition of / 
beyond a, contradicting the maximality property of the interval (a, ft). 

Suppose there are sequences (tj)jen, ( s j)jeN C (a, /3) such that tj — > a and 
Sj — > a and /(tj) — > a while /(sj) — > b as j — > 00. Then for any rj' e (a, 6) there is 
a sequence (rj)jeN C (a, /?) such that rj — > a and /(rj) — > 77', whence 77' S a(H a ). 
This would imply that (a, 6) C a(H a ), which is impossible. □ 

We next turn our attention to the question of Lipschitz-continuity of the func- 
tions fj in Lemma 2.1. With §v ■ [0,1] — > [0, 00) as in (A. 3), we study potentials 
from the classes 

V a := {V e V I 3C > 0: «V(s) < Cs Q ,V0 < s < 1}, (A.5) 

where < a < 1. The class P a consists of all periodic functions V € "P which 
are (locally) a-H61der-continuous in the Li-mean; for a = 1 this is a Lipschitz- 
condition in the Li-mean. The class "Pi is of particular practical importance since 
it contains the periodic step functions. It will be shown below that V\ coincides 
with the class of periodic functions on the real line which are locally of bounded 
variation. We first prove Lipschitz-continuity of the eigenvalues of H t for V £ V\ . 

A. 3. Proposition. For V € V\, let (a,b) denote any of the gaps Lfe of H and 
let fj: (aj,/3j) — > (a, b) be as in Lemma 2.1. Then the functions fj are uniformly 
Lip schitz- continuous. More precisely, for each gap T k there exists a constant Ck > 
such that for all j 

\fj(t) - fj{t')\ < C k \t - t'|, a 3 <t,t'< fr. 

Proof. As in the proof of Lemma 2 .6 we can find a finite number of levels E\ , . . . , Et <G 
(a, b) and a partition = t < n < . . . < T£_i < 17 = 1 such that Ej £ cr(H t ) for all 
t € Ij := [Tj-i,Tj] and for j = 1, ...,£. Now V e Vi implies \W t - W^'lli.iocunif = 
i?y(t — t') < C\t — t'\ and we conclude with the aid of Lemma A. 2 that there are 
constants c\ , . . . , ct > such that 



\\(H t - Ej)- 1 - (H r - Ej)- 1 ]] < Cj \t-t'\, t,t' & Ij. 
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This implies that the min- max- values Hk(s) of (H s — Ej) 1 satisfy 

IMt) -Vk{t')\ <cj\t-t'\, t,t' eij. 

By the spectral mapping theorem, the eigenvalues of H t in (Ej, b) are in bijection 
with the eigenvalues of (H t — Ej) -1 in { b } E . , oo). We now let C :— max{q, . . . , c e } 
to finish our proof. □ 



A. 4. Remarks. 

(a) By the same argument, we obtain the following result on Holder-continuity: 
If < a < 1 and V € V ai then each of the functions f 3 ■ : (ctj,(3j) — > (a, b) 
is locally uniformly Holder-continuous (as defined in [GT]), i.e., for any com- 
pact subset C (aj,(3j) there is a constant C = C(j,a'j, (3'j) such that 
\fj(t) - fj{t )\ < C\t - t'\ a , for all t,t' e [a'^fy]. Note that our method docs 
not necessarily yield a uniform constant for the whole interval (otj ,(3j), much less 
a constant that would be uniform for all j. 

(b) For analytic potentials V, it is shown in [Kl] that the eigenvalue branches 
fj in Lemma 2.1 depend analytically on t. This is a simple consequence of the 
fact that, for real analytic V, the Ht form a holomorphic family of self-adjoint 
operators in the sense of Kato. In [K2], the author proves that the fj are squares 
of W\- functions near the gap edges if the potential is in £2(1). 

We finally give a characterization of the class V\ . 

A. 5. Proposition. Let BV\ oc (M.) denote the space of real-valued functions which 
are of bounded variation over any compact subset of the real line. 
Then Vi =VDBVi oc (R). 

It is easy to see that any 7ePfl BV\ oc (R) belongs to V\: certainly, any VeP 
which is monotonic over [0, 1] is an element of V\ and any function of bounded 
variation can be written as the difference of two monotonic functions. 

The converse direction is established by the following result due to J. Voigt, 
Dresden; cf. also [EG, Chapter 5] for related material on BV^-functions of several 
variables. 

A. 6. Lemma. Let f G L\\ oc (M, R) be periodic with period 1 and suppose that 
there are c > 0, e > such that 

[ \f(x + t) - f(x)\dx < ct, W0<t<e. (A6) 
Jo 

Consider f as a function in Li(T), with T denoting the one- dimensional torus. 

We then have: the distributional derivative df is a (signed) Borel-measure [i on 
T and there is a number a G R such that f{x) = a + fi([0,x)), a.e. in [0, 1) ~ T. 
Ln particular, f has a left- continuous representative of bounded variation. 
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Proof. Defining r\ : C 1 (T) -> R by 

77(^5) := - / y'/d.r. 

we may compute 







- / <p'f dx = lim / -t)- <p(x))f(x)dx 

Jo *^°Jo 1 

f 1 1 

= |imy ^)-(/(.x + t)-/(x))dx, 

and the assumption (A. 6) yields the estimate |»7(v?)| < c IMIco- Since C 1 (T) is 
dense in C(T), the functional 77 has a unique continuous extension to all of C(T); 
we denote the extension by the same symbol r\. By the Ricsz representation theorem 
there is a measure \i such that rj(ip) = J <pdfi for all ip e C(T). Furthermore, for 
ip € C 1 (T) wc have — tp'fdx = pdfi, and we see that /j = 9/ on T in the 
distributional sense. The choice ip := 1 yields J T d/j = —f^tp'/dx = and the 
function f(x) := /x([0,x)) satishes 9/ = /i. This is easy to check: for ip e C 1 (T) we 
have 



f(p'dx= / d/^(y)93'(x)dx 

J JO Jo<y<x 

= I / <^'(a;)da;d/i(y) = - / 

J0<y<lJy J[0, 



0<y<x 

ip(y) dfj,(y). 



M) 

We therefore see that d(f — /) = 0; hence there is some a such that / — / = a. □ 
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